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Abstract: Beginning with a special form of the Einstein- Rosen metric, 
we find new cosmological solutions of the Einstein equations, having two 
hypersurface-orthogonal Killing vectors , with ideal fluid. The equation of 
state is in the most cases of the form p = jp. 



1 Introduction 



The actual universe is neither exactly spatially homogeneous nor exactly 
isotropic. However, it is much too complicated to implement these real fea- 
tures of the universe into our models. The cosmological principle is the sim- 
plifying assumption which leads us to a very simple model having only one 
degree of freedom: the cosmological scale parameter of FRW-metrics. Now, 
the shortcommings of the standard cosmological model is well known, and 
even the old-, new-, choatic-, and extended- inflationary models are not free 
from problems. A better understanding of the standard model, and even the 
inflationary ones, towards a more realistic model of the universe is achieved if 
one could study solutions without the implementation of cosmological prin- 
ciple. 

Anisotropic spatially homogeneous cosmological models have been stud- 
ied extensively [1]. Even questions of inflation and bubble nucleation in 
anisotropic models have been considered [2]. However, the inhomogeneity 
is not as easy as anisotropy to implement. This is why we do not know much 
about it yet. Considering the most general inhomogeneous model is yet out 
of question, and one should therefore try to understand the simplest inho- 
mogeneous models. These are models with one-dimensional inhomogeneity 
having two space-like commuting Killing vectors known as orthogonally tran- 
sitive G2 cosmologies [3-5]. Up to now, there are few solutions of this type 
known which shows a variety of features [6- 10]. The first class of solutions 
was given by Wainwright and Goode[6]. The class of solutions found by Ruiz 
and Senovilla[8] have quite different behaviours as far as the structure of sin- 
gularities is concerned. There are solutions with a big-bang-like singularity, 
solutions with time-like singularities only in the Weyl-tensor, solutions with 
singularities in both the Ricci- and the Weyl-tensor, and also singularity-free 
solutions. Aguirregabia et.al. [9] are interested in inflationary behaviour of 
inhomogeneous models with exponential-potential scalar field as the source, 
and they have found solutions with multiple inflationary behaviour for which 
the deceleration and the inflationary phases interchange with each other sev- 
eral times during the history of the model. Although these models tend to 
homogenize, but they never isotropize. 
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2 Metric and field equations 



The simplest one-dimensional inhomogeneity can be formulated via the Einstein- 
Rosen metric[4,5]: 

ds 2 = -Adt 2 + Bdx 2 + C{Ddy 2 + D~ l dz 2 ) (1) 

where A, B, C, and D are functions of t and x. This metric represents a space- 
time having two commuting hypersurface -orthogonal Killing vectors. To 
simplify the solution all the authors assume the unknown functions to be 
separable in t and x. Besides this, authors of [6] and [9] assume D = D(t) 
and C = C(t), respectively. Assuming A = B, and that the time part of 
A be a power of that of the function G we obtain, in essence, the class of 
solution of [8]. 

We have been able to find some classes of new solutions of the metric in 
the following form 

ds 2 = -F 2 dt 2 + T 2m dx 2 + T.G(T n Wdy 2 + T~ n W~ l dz 2 ) (2) 

Here is T a function of time and F,G, and W functions of x only. We assume 
the universe to be filled with perfect fluid, having an equation of state to be 
determined later on. The Killing vectors are and J^; the unit velocity 
vector of the fluid is 

.— (3) 

The energy momentum tensor of perfect fluid is taken to be 

Tp, = (p + p)u^u u + pr]^ (4) 
To calculate the field equations we use the orthonormal tetrad 

J = Fdt , u x = T m dx 



v = T (n+1)/2 VGWdy , u z = T^- n)/2 ^G/Wdz (5) 

Now, for the line element (2) and the matter source (4) we obtain the fol- 
lowing non-trivial and independent Einstein field equations: 
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The equation (6) can be integrated to lead 

F 2 = Ql-2m w n 



= (8) 

(9) 
up (10) 

(11) 



In equation (7) the variables can be separated. We call the separation pa- 
rameter eb 2 , where e is choosen to be or ±1: 
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Using equation (11) and defining new variables 



W n G' 
a= W ' P=G 

the field equations (6-8) can be written in the following form 

F' n 1 — 2m „ 
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F 2 2 1 



(12) 
(13) 
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The four equations (14-17) can be solved to find the unknown functions F, 
T, G, and W. The equations (9) and (10) then give pressure and density 
as functions of x and t. It remains to be verified if there exists an adequate 
equation of state. 

The kinematical quantities can now be calculated easily. The rotation 
comes out to be zero. The other non-zero components of expansion, acceler- 
ation, and shear, in tetrad components, are 
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3 Solutions 

We differenciate the cases according to the different values of e. 
Case I: e = 0. 

Equation (17) can then be integrated immediately to give T: 

T(t) = C 2 [C 1 + (m+l)t] m+1 
where C\ and Ci are constants of integration. Now assuming 

G = W k 



(23) 



(24) 
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with constant k, we see that (3 = ha, where /3 and a are defined in (13). 
With this ansatz the equations (15) and (16) can be solved to give 



n ± \/n 2 + 3 - 4m 2 o ^ 

fc = ^ , n - 4m + 3 > (25) 

4m - 3 

W = C 2 [C[ + (n + fc(3 - 2m))x]™^) ( 2 6) 

From (24) and (11) the other metric functions are determined. The remaining 
field equations (9) and (10) gives p and p and the equation of state. It comes 
out 

Kp = Kp = *0Z*±± W -W~* (^—) 2 (27) 

Hence, the fluid obey a stiff equation of state. For the density to be positive 
there must be 

4m - n 2 + 1 > (28) 

This gives together with (25) 

4 > 4m - n 2 + 1 > (29) 

The case m = — 1, which would lead to constant T, is now automatically 
excluded. Obviously the zeros of the functions T(t) and W(x) leads to sin- 
gularities. The former is a big-bang-like singularity. The latter one, which 
is generic in most of inhomogenous solutions yet found has not been studied 
throughly. To study such type of singularities, which we call them 'wall-like 
singularity', it is necessary to pick up some tractable special solutions suit- 
able for physical interpretations. 

Case II: e = +1. 

The integration of (17) in the following cases is straightforward: 

m = : T{t) = C x e u + C 2 e- bt (30) 



m = +1 : T(t) = ±\Jb 2 t 2 + dt + C 2 (31) 
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Now, assuming 

T(t) = [U(t)\* (33) 
with constant A, the equation (17) can be integrated in general to give 

If 2 = C 1 + b 2 {m+l) 2 U^ (34) 

For the integration to be possible we have choosen A to be equal ^zj- Ci 
is again a constant of integration. For C\ — 0, the equation (34) can be 
integrated. We obtain 

T(t) = (C±mbt)™ with m^0,-l (35) 

Note however, that for m = +1 this solution is a special case of (31). Now 
we find two classes of solutions for the x-dependent functions. 

Case Ila: Let 

W = Ae qx (36) 
G = Be sx (37) 
where A, B, q and s are constant. From (14-16) we then find 

F = ^M 1 -^/ 2 (38) 
^ e[nq+ (i-2m)s] x j^ g 2 + ^^qs^j = nb 2 (39) 

2 

+ 2m(2n - 9 1) - n + 3 g S ] = -b 2 (2m -n-l) (40) 

From (39) and (40) we must have 

nq+{l-2m)s = (41) 

The above equations finally lead to 

^ 2 2?n 

-g 2 + — ^— gs = n6 2 /^ 2 (42) 
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-n 2 + n - 2 9 , . , 2m(2n - 1) - n + 3 
g 2 + m(l - 2m)s 2 + ^ ^ qs = 

= -(2m - n - l)b 2 /(i 2 (43) 

The parameters q, s and are obtained from (41-43). For the unknown 
functions we get 

F = fjL , n = A n/2 B^ 2m ^ 2 (44) 

W = Ae ± ^ V ^ =Tx (45) 

G = Be ± Tv^f x (46) 

As we are working with real metric, the above solutions are restricted to 
m > 1/2. Hence in calculating p and p we have just the choise of (31) and 
(35) for T(t). Taking the general choice (35), we obtain from (9,10): 

b 2 (m+ l)(2m-n 2 - 1) 1 
Kp ~ V 2m~^~l (C 2 ±mbty { ' 

V (rn-l)(2m-n 2 -l) 1 

y 2/i 2 2m -1 (C 2 ±mbt) 2 1 ; 

with m > 1 and n 2 < 2m — 1. The equation of state comes out to be 

Tfl — 1 

P = IP , 7 = —r , < 7 < 1 (49) 

m + 1 

Taking the choice (31) for m — 1, we obtain 

6 2 

«p = — (1 - n 2 )x 
P 



x 



^ + Cl(+ (iz" 2 )^-W + i )c ^ 



165 2 (l — n 2 ) j + + (*)"* (50) 

KP = ( Cl2 " 462C2 ) ( 62t2 + Cl * + ^ (51) 

The proportionality factor between pressure and density is now a function of 
time. Looking for a simple equation of state of the form 

P = IP ( 52 ) 
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with constant 7, we are led to the following conditions 

d = 2bC , C 2 = C 2 
This leads to the following equation of state: 

b\l-n 2 ) 1 

p = ' Kp = f (bt + c-Y (53) 

This correspond to the solutions (47, 48) for m = 1. Hence, we have a dust 
solution with a standard big-bang-like singularity. 

Case lib: Let 

a =m • p= w> (54) 

with constant D and E. Substituting in the equations (14-16) we obtain 

F(x) =k + A x (55) 

W(x) = h(k + Aqx) d ' Ao (56) 
G(x) = k 2 {k + A x) E / A ° (57) 

where 

Ao = ±nD + ±(3-2mE)-^ (58) 
The following consistency relations between constants must also hold: 



k 2 = k^- 1 (59) 

ED = nb 2 (60) 
E 

- D 2 + (3 - Am)E 2 + InDE + 2n(m - 1)6 2 — + (2m - n 2 - l)b 2 = (61) 
There are two solutions of (60, 61): 

D 2 = n 2 b 2 (62) 

D 2 = (2m - l)b 2 (63) 
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Choosing (63) we obtain from (58) that Aq = 0, which is unacceptable due 
to the relations (56) and (57). Therefore we try (62). Using (35) for T(t), it 
comes out 

Kp — Kp — (64) 

Now we have still the choices (30-32) at our disposal. The case m = leads 
to the following relations for pressure and density: 



Kp = Kp 



[n 



l)b 2 



■C 1 C 2 (C 1 e bt + C 2 e- u )' 



(k + A x) 2 

The condition for the positivity of the density is 

{n 2 - \)C X C 2 > 
The case m — — 1 leads to negative density. But for m — 1 we get 



(65) 
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Kp = Kp = 

1 (b 2 t + Cx/2) 2 



b 2 b 2 t 2 + dt + C 2 



1 



(k + A x) 2 (b 2 t 2 + Cit + C 2 ) 



(67) 



The equation of state is again in the stiff form. 
Case III: e = — 1. 

Solutions for this case can be find using again the forms 

W = Ae qx 



and 

This choice leads to 



G = Be s 



W = Ae^^ 

_ rtb 1 
Q = Be T * VT^2m 

F = p , p = A n/2 S (1_2m)/2 , m < 1/2 
The equation (17) leads to 

T(t) = (C ± imbt)™ m^O 



(68) 

(69) 
(70) 

(71) 
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For m = one gets a sinosoidal function which leads to violation of energy 
condition and is therefore unphysical. We obtain real values of (71) for 

m = 47 

this leads to 



lez-o 



(72) 



T(t) = (^) 4 < 

Pressure and density are calculated from (9) and (10): 
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21 - 1 

with the equation of state in the form p = •yp. The positivity of energy leads 
to / > 0; but the energy condition < p < p leaves no choise for n and 
/. Seeking another solution along the lines of the case lib leads to complex 
metric functions. 



4 Discussion 

The solutions we have found shows very diverse singularity behaviours. There 
are some with big-bang- and wall-like singularities, such as (27) and (67). 
Others show just big-bang-like singularities, such as (47,48), (50,51,53), and 
(74,75). The class of solutions leading to (65) is peculiar. It shows no big- 
bang singularity. In fact, there are some solutions with no singularity in time. 
However, it has a wall-like singularity. The nature of all these singularities 
should be studied very carefully. One has to consider special solutions and 
look for any physical interpretations of them. We are considering some in- 
teresting solutions of this class and those found by other authors, which will 
be subject of other publications. 
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